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Abstract. 

Ultra-high energy cosmic rays (UHECR) above an energy 
threshold of tens of EeV might undergo only small deflec- 
tions due to interstellar magnetic fields. Their arrival direc- 
tions would then point to regions of possible hadronic accel- 
eration processes, which are likely to be also sources of high 
energy neutrinos. To search for such cosmic accelerators, 
we present here the first high statistics analysis of directional 
correlations between neutrino candidates from the IceCube 
Observatory and UHECR events. Data taken with IceCube 
in a configuration of 22 strings provided the basis for using 
published events from both the Pierre Auger Observatory and 
the HiRes experiment as reference directions in a search for 
coincidences with neutrinos. The analysis was optimized ac- 
cording to strict blindness criteria and showed an excess of 
neutrino candidates close to UHECR directions with a proba- 
bility of 1 % to occur as a random fluctuation, consistent with 
a background-only hypothesis. An extension of this analysis 
to include newer IceCube data, taken with 40 strings and us- 
ing a likelihood analysis, is discussed in the outlook. 



1 Introduction 

The successful operation of the HiRes experiment until 2004 
and the growing data sample provided by the Pierre Auger 
Observatory (PAO) made high statistics studies of UHECRs 
possible. Nevertheless, their sources remain unknown and 
the composition features have not yet been conclusively de- 
termined ("AbrahanTetig |20T0) |Ab^^ The 
limited knowledge about astrophysical magnetic fields pro- 
hibits detailed predictions about their propagation, but at en- 
ergies above 10^^ eV, at least protons and light nuclei should 
undergo only small deflections on the order of a few degrees. 
Consequently, these cosmic rays will point back to regions 



of the sky where astrophysical sources achieve hadronic ac- 
celeration up to the highest energies, and are thus promising 
places to search for high energy neutrino emissions. 

Up to now, no cosmic point sources of neutrinos at TeV 
or higher energies have been identified. Searches for these 
neutrinos are performed on data samples dominated by atmo- 
spheric background events and in consequence their discov- 
ery potential can be boosted by limiting the trials via testing 
only directions defined by UHECR events. Such a multi- 
messenger approach has the primary objective to search for 
efficient neutrino emission connected with hadronic acceler- 
ators but could ultimately help to learn more about the distri- 
bution of UHECR sources and astrophysical magnetic fields. 



2 Instruments and data samples 



The IceCube Observatory (Karg 20I0| l is currently the 
largest neutrino telescope and was used in this UHECR cor- 
relation search. Located in the deep ice at the geographic 
south pole, photomultipliers in digital optical modules ( ,Ab^ 
basi et al. 2009d[ l are arranged on strings, carrying 60 mod- 
ules each and reaching down to 2450 m below the surface. 
The array serves to measure Cherenkov light from charged 
particles, in particular secondaries from neutrino interac- 
tions, and a nanosecond resolution of photon arrival times 
allows reconstructing the direction of muon tracks. The de- 
ployment of a total of 86 strings in a volume of a cubic kilo- 
meter will be finalized until the beginning of 20 11. 

The first UHECR correlation study has been performed 
with IceCube data collected over 276 days in 2007/2008 with 
a configuration of 22 strings. It was based on the muon track 
sample selected for the first point source search that included 
not only upward going muon tracks, but also events from 
above the horizon, a region previously excluded due to in- 
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creased background fluxes ( ,Abbasi et al. 2009c I. The events 
from the northern sky are mainly induced by atmospheric 
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Fig. 1. Arrival directions of UHECR in equatorial coordinates, with 
the galactic plane shown as a curved line. 13 events are from HiRes 
(blue), 27 from the PAO 2007 sample (red) and 42 additional events 
(orange) released in 2010. 



neutrinos, while the southern sky is dominated by downward 
going atmospheric air shower muons. The neutrino energy 
range covered with this selection for a hypothetical sig- 
nal spectrum reaches from TeV to PeV for upward going 
events, limited by absorption in the Earth. The energy range 
for downward going tracks covers PeV to EeV, because lower 
energies are suppressed through background reduction cuts. 

For these correlation searches, publicly available UHECR 
data from two extensive air shower detectors are used as 
reference directions. One is the HiRes experiment, operat- 
ing between 1999 and 2004 in stereo mode with two fluo- 
rescence detectors (Abbas i et al.| |2009a| l. The other is the 
southern PAO array which is taking data since 2004 with a 
hybrid technique of fluorescence telescopes and Cherenkov 
detectors. The first data on arrival directions of UHECRs 



published by PAO in (Abraham et al. 2008^ comprised 27 
events, collected between January 2004 and August 2007. 
The chosen energy threshold of 56 EeV was optimized in a 
search for correlations with nearby AGN together with two 
other search parameters. In that analysis, the hypothesis of an 
isotropic distribution of UHECRs was rejected with a prob- 
abihty of 99% and a correlation with the AGN sample was 
claimed. Recently, the results were updated in ( |Abreu et al.[ 
20 10] ! with 42 additional events detected through the end of 
2009. The degree of correlation with the AGN catalog has 
decreased for the full sample, but this has no impact on the 
neutrino analysis presented here, since only the reported co- 
ordinates of UHECR arrival directions themselves are con- 
sidered. Following the PAO reports, the HiRes collaboration 
applied the same energy threshold as in the PAO study to all 
their events collected in stereo mode. They found 13 events 
and published the arrival directions in (Abbasi et al. 2008[ l. 
HiRes could neither confirm a correlation with the AGN cat- 
alog nor find an indication of anisotropy in a separate test 
with a different event selection. 

At the time of our first correlation search, only the orig- 
inal 27 PAO and 13 HiRes events had been made public. 
Because of the low signal expectations, we did not use Ice- 



Cube 22 strings data with declinations below —50°, leaving 
35 UHECR arrival directions. The new IceCube event se- 
lection for 40 strings covers all declinations and thus can be 
correlated with the full sample of 82 UHECR events shown 
in Fig. [1] 

3 Binned correlation method 

A possible correlation between UHECR directions and the 
IceCube data can be tested by counting neutrino candidates 
in circular bins of a pre-defined radius centred on the arrival 
directions of the charged cosmic rays. In ( |Petrovic | |2009[ ), 
an application of such an approach to generic neutrino and 
UHECR data samples is discussed. To calculate the signif- 
icance of the observation, the sum of events over all bins 
has to be compared to the expectation from pure background 
samples. These can be produced by randomizing the indi- 
vidual right ascension values of measured data, because the 
acceptance varies only with declination due to IceCube's lo- 
cation and the integrated live time of a year. Counting Ice- 
Cube events over all bins is performed by testing if the min- 
imum angular distance to any UHECR direction is smaller 
than the fixed bin radius to avoid double counting in overlap- 
ping bins. The distribution of results from scrambled samples 
is then integrated over all event numbers larger than the one 
observed. This yields the p-value, defined as the probability 
for the measured outcome to appear in a background-only 
scenario. 

Simulated sky maps with added signals are required to as- 
sess the discovery potential and optimize the analysis for re- 
alistic expectations. Since only the UHECR directions can 
be used as reference positions for simulated neutrino sources, 
the effects of magnetic fields on the ionized nuclei have to be 
taken into account. There are few established facts about as- 
trophysical magnetic fields and predictions vary concerning 
their strengths and orientations. The PAO collaboration as- 



sumes deviations around 3° and up to 7° in ([Abraham et al. 



2008j l and references therein. Other authors like ( Kachel-| 
riess et al. 2006| l also predict average deflections to be on 
the order of a few degrees in the energy range under con- 
sideration, though possibly stronger deviations occur close 
to the galactic center Extensive simulations for galaxy clus- 
ters discussed in ( |Dolag et al.j [2005) ) show that the impact 
of extra-galactic fields on UHECR propagation is negligible. 
These predictions are only valid for proton UHECRs, and 
deviations for heavier nuclei are expected to be significantly 
larger which would then strongly limit the prospects of a neu- 
trino correlation approach. 

To account for unknown magnetic deflections in the sim- 
ulation, neutrino sources were individually shifted with re- 
spect to the UHECR directions. Based on the above ref- 
erences, a width of 3° was chosen for a Gaussian function 
according to which randomly orientated angular displace- 
ment values were generated in each trial individually for all 
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UHECR directions. The angular resolution of the air shower 
reconstruction for the chosen data is better than 0.9° ( |Abbasi| 
|et al.| |2008[ |Abreu et al.| |2010| l and was thus considered as 
negligible compared to the deflection uncertainty. 

In the absence of any information on individual muon neu- 
trino expectations for the hypothetical sources, equal point 
source fluxes from all UHECR directions, following an E^^ 
spectrum, were assumed in the simulation. Accounting 
for the variations of IceCube sensitivity with declination as 
shown in ( Abbasi et al.[ |2009c) , average numbers of neutri- 
nos were derived for each source direction and used to gen- 
erate event numbers in each trial according to Poisson distri- 
butions. Simulated signals were placed into scrambled sky 
maps by distributing them around the shifted source posi- 
tions with angular deviations that followed the IceCube point 
spread function. 

The distributions of binned event counts for simulations 
with and without signal can be well fitted with Gaussian 
functions. They allow deriving a 5a (one-sided Gaussian 
tail) threshold from the background scenario and then deter- 
mining the flux strength for which the probability of such an 
excess reaches 50%, called the discovery potential. The bin 
radius is the only free parameter in this correlation approach. 
For 22 strings data, various radii were tested with respect to 
the best discovery potential for simulated signals and an op- 
timal value of 3° was found. Simulations with an increased 
width for the source shifts showed that the optimal value does 
not become significantly larger even for greater deflections 
due to the additional background covered by the bins. De- 
tails on the optimization are provided in ("Lauer 2010"). Since 



the prescription was fixed before applying it to the true event 
coordinates, no trial factor enters into the final p-value. 

The resulting discovery potential for the 22 strings search, 
expressed as an E^'^ muon neutrino flux normalization for 
d^/dE = <^>oE-'^, is $0 = 1.0- 10"^ GeV cm^^ g-i. This is 
approximately a factor of three smaller than the upper limit 
on the average muon neutrino point source flux for 35 equal 
sources derived from the best diffuse neutrino Umit in the 
same energy range in ( [Ackermann et aL)|2008] ). 



4 Results from the search on IceCube 22 strings data 

After applying the analysis prescription to the IceCube data 
by scanning the 35 bins of 3° radius, a total of 60 events 
were observed. The mean background expectation based 
on scrambled sky maps was 43.7 events. The probabil- 
ity for an excess of the same or a larger magnitude within 
the background-only scenario is 0.0098, equivalent to 2.3(t. 
This result is compatible with a background fluctuation. No 
particularly strong clustering of IceCube events on smaller 
scales than the bin size was observed and the muon energy 
estimates for the 60 events featured no unusual deviations 
from the full distribution of background data. 



Under the benchmark emission hypothesis of equal neu- 
trino fluxes from all 35 R^O and HiRes UHECR direc- 
tions, the result can be used to derive an upper limit on 
the muon neutrino flux normalization per source, yielding 
"J'lim — 0.9 • 10^^ GeV cm^^ s^^ for a mean magnetic de- 
flection of 3°. Since this uniform emission scenario is only 
a convenient framework for interpreting the event counts, 
the limit should be considered at best as an average level of 
probed flux strengths. 



5 Outlook for the search with IceCube 40 strings data 



The newly released PAO arrival directions from ( Abreu et al. 



20101 allow us to extend the analysis to a total of 82 UHECR 
reference directions. Based on this a search on the indepen- 
dent IceCube data collected with 40 strings over 376 days is 
in progress. The event selection is adopted from the all-sky 
point source search as presented in ( [Abbasi et al.j |2010b[ ), 
comprising 36,900 muon tracks consistent with expecta- 
tions from atmospheric background simulations. Scrambled 
sky maps and simulations with added neutrino signals from 
shifted sources were produced in the way described above to 
study the performance of the binned search as well as a new 
unbinned likelihood approach. 

The differences in the point source event selection for 40 
strings compared to the 22 strings sample made it neces- 
sary to apply an additional cut on the energy estimator to 
perform the binned search in a similar way. The cut was 
optimized for sensitivity and is parametrized as a zenith de- 
pendent threshold, decreasing from the downward to the up- 
ward going regime. It resulted in 6013 remaining neutrino 
candidate events. Applying the binned search method to the 
simulations proved that the change of the discovery poten- 
tial with bin radius is again small. Hence it was decided to 
keep the value of 3° for consistency with 22 strings. The 
resulting discovery flux threshold for the binned search is 
$iim = 0.7 • 10^^ GeV cm-2 s'^ The impact of the addi- 
tional UHECR directions is small, because most of them lie 
in the southern hemisphere where the E^^ acceptance is lim- 
ited due to the high energy threshold. 

To improve the sensitivity by including individual event 
parameters in the probability calculation, we adopted the un- 



binned likelihood method from (Braun et al.| |20()8| to the 
correlation search. Its use in IceCube point source searches 
has been established in ( [Abbasi et al.| |2009b| l. The like- 
lihood formula Eq. ([T} is a product over the contributions 
with index i from all IceCube events. As done in ( ,Abbasi| 
et al. 20 10b I for stacking searches on astronomical object 



classes, we make the signal term a sum over the contribu- 
tions from many reference points, in our case the 82 UHECR 
arrival directions Uj. The expectation for the average num- 
ber of neutrinos which would be detected from direction j for 
a common source spectrum E ' is included via the accep- 
tance weighting function R{j,'-f). It depends on the position 
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of each hypothetical source and corrects for the declination- 
dependent energy response of the detector The sum over the 
UHECR events with index k serves as normalization. 

ntot , 

I J k 

+ (i-J^\b{^.,.e,)) . (1) 

V "tot / 

By minimizing the negative logarithm of C, the best fit val- 
ues for the total number of signal events Us and a common 
spectral index 7 is determined. The comparison of the like- 
lihood value with that for Ug = yields the p-value for the 
null hypothesis. The background probability density func- 
tion (PDF) B{xi,Ei) depends on the individual event energy 
estimator Ei and the declination of its location Xi and is de- 
rived directly from data. For the PDF Sj{xi,(Ti,Ei,j) repre- 
senting the signal hypothesis with respect to direction j, the 
energy dependence is obtained from neutrino simulations for 
spectral indices 7 between 1 and 4. In the spatial part, the 
point spread function is represented by a Gaussian function 
with a width ai that is obtained individually as a direct re- 
construction error estimate ( |Abbasi et al.||2010 F). This func- 
tion has to be convoluted with the probability distribution for 
a source shift due to magnetic deflection, for which we use 
again a Gaussian with a width of fr^ = 3°, as discussed with 
respect to the simulation above. This convolution results in a 
broadened Gaussian function for the spatial PDF, 

(2) 

The algorithm is applied to sky maps of the full 40 strings 
sample, with simulated neutrinos injected from shifted 
sources as outlined above. The resulting discovery potential 
for the scenario of equal E~'^ fluxes from all 82 directions 
is (f)iijn = 0.2 • 10^^ GeV cm^^ g-i xjjis is equal to an av- 
erage of 52 signal neutrinos in the point source sample of 
36,900 events. The improvement over the binned method is 
largely achieved through the direct use of event energy es- 
timators and individual reconstruction errors and the corre- 
sponding PDFs. The acceptance function R gives additional 
discrimination power by weighting dkections according to 
their hypothetical contributions, but this will only be a ma- 
jor advantage if the assumption of many more or less equal 
sources is justified. 

Testing the same algorithm on simulations of neutrino 
sources with the width of the source shift distribution in- 
creased to 6° raises the discovery flux threshold by approx- 
imately 50%. This reflects a good degree of flexibility, in 
particular when compared to the binned search where the in- 
creased shifts worsen the discovery potential by a factor of 
2. 



6 Conclusions 

We performed the first large statistics analysis of directional 
correlations between 35 UHECR arrival directions reported 
by PAO and HiRes and IceCube data collected over one year 
with 22 strings. The observed excess has a 1% chance to oc- 
cur as a background fluctuation and is consistent with the 
background-only hypothesis. The IceCube data collected 
with 40 strings and newly released UHECR directions are 
currently being used to perform an extended search, employ- 
ing a likelihood approach to improve the sensitivity com- 
pared to the binned event counting. 
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